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A commentary on 

Multilevel control of run orientation in 
Drosophila larval chemotaxis 

by Gomez-Marin A., and Louis M. 
(2014). Front. Behav. Neurosci. 8:38. doi: 
10.3389/fnbeh.2014.00038 

In their article, Gomez-Marin and Louis 
(2014) found that runs in Drosophila lar- 
val chemotaxis bend toward the direc- 
tion of higher concentration. This steering 
process called weathervaning or klinotaxis 
was previously discovered in the worm 
C. elegans (Ward, 1973; lino and Yoshida, 
2009) and raises the question: how is it 
performed? Although Gomez-Marin and 
Louis (2014) found that run orientation 
relies on the detection of the lateral gra- 
dient component, the underlying con- 
trol mechanism is largely unknown. Two 
alternative strategies might be employed, 
direct error correction and proportional 
navigation (Figure 1A). In direct error 
correction, the animal would turn pro- 
portionally to the error; that is, the local 
bearing angle between the current head- 
ing and the desired direction to the goal 
(intensity peak). In proportional naviga- 
tion (Murtaugh and Criel, 1966; Zarchan, 
2012), the objective is to maintain a con- 
stant line of sight (LOS) angle during 
motion (Figure 1A). This tactic has been 
known by sailors for many years as a mean 
to detect collision courses and is currently 
implemented in guided missiles. Thus, 
in proportional navigation, the animal 
would turn proportionally to the deriva- 
tive of the LOS angle. Here we argue that 
klinotaxis, as observed in Gomez-Marin 
and Louis (2014), does not result from 
direct error correction. Instead, klinotaxis 



is directed through proportional naviga- 
tion that is a basic form of navigation 
(where the goal direction is estimated 
using gradient sensing rather than spatial 
memory) involving indirect error correc- 
tion through an attempt to keep the LOS 
angle constant. Our claim is supported by 
theoretical analyses as well as computer 
simulations. 

WHAT KIND OF MOTOR OUTPUT 
RESULTS FROM DIRECT ERROR 
CORRECTION AND PROPORTIONAL 
NAVIGATION? 

The agent, either Drosophila larva or C. ele- 
gans, is modeled as a single point (x, y) 
moving at constant speed v. The equations 
of motion are given by 



dx 
dt 
dy 
dt 



vcos y 
vsin y 



(1) 



where y is the direction of movement with 
respect to a global reference (Figure 1A). 
The agent aims at orienting toward the 
goal by controlling its turning rate to = 
dy/dt. The orientation error is given by 
the local bearing angle 8 = a — y where 
a is the line of sight (LOS) angle; that 
is, the absolute direction of the bee- 
line from the agent to the goal. What 
kind of motor output may result from 
direct error correction and proportional 
navigation? 

A direct error correction strategy would 
produce a motor output (turning rate to) 
proportional to the error (local bearing 
angle B) 



with proportional constant K. Yet, no 
evidence of such a linear relationship 
was found in (Gomez-Marin and Louis, 
2014). Instead, Figure IB clearly shows 
that co ~ KsinB. A similar turning rate 
function is observed in C. elegans klino- 
taxis (see Figure 2A in lino and Yoshida, 
2009) and in other types of taxis, e.g., 
gyrotaxis (Codling et al., 2008). Can such 
sinusoidal reorientations be obtained with 
proportional navigation? Proportional 
navigation is a dedicated term for a 
basic form of navigation that aims at 
maintaining the LOS angle a constant dur- 
ing motion (Murtaugh and Criel, 1966; 
Zarchan, 2012). It was originally devel- 
oped for missile guidance toward moving 
targets. Recently, it has been employed for 
navigating a robot toward a fixed goal with 
the kinematics given by Equation ( 1 ) and 
was defined as follows (Belkhouche and 
Belkhouche, 2007). 



K- 



r da 
dt 



(3) 



so that the agent moves in a straight 
line when a remains constant. In the 
Supplementary Materials, we show that 
the rate of change of a is given by v sin 8/r 
so that 



co = K- sin/3 
r 



(4) 



co = KB 



(2) 



Thus, proportional navigation is expected 
to produce a motor output co proportional 
to the sine of 8, similar to that observed 
experimentally (Figure IB). Note however 
that Equation (4) is useless in practice as 
both the distance r to the goal and local 
bearing angle 8 are unknown to the agent. 
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Error correction vs proportional navigation 

Turning rate: 
dy 

co — — — 
dt 

Error correction: 
co = K.p 




Goal 

(intensity peak) 



Proportional navigation: 

da 
dt 



co = K^ 
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FIGURE 1 | (A) Error correction vs. proportional navigation. The turning 
rate is proportional to the local bearing angle in error correction and to the 
rate of change of the LOS angle in proportional navigation. See text for 
the details. (B) Dependency of the turning rate on the local bearing angle. 
The experimental data (in blue) are reprinted with modification from 
Figure 1C in Gomez-Marin and Louis (2014). The straight line represents 
error correction (Equation 2 with K = 1.5). The dashed curve is for 
a> = Ksinfi with K = 1.5. (C) Longitudinal sampling of the gradient along 
the direction of motion while the animal moves forward. See text for the 
details. (D) Lateral sampling of the gradient perpendicular to the direction 
of motion during head sweeps. See text for the details. (E) Dependency 
of the turning rate on the local bearing angle for simulated proportional 
navigation (black curve, mean ± s.d., n= 100 trajectories) as compared to 
experimental data (blue curve, same as in panel B). The inset shows a 



superposition of the trajectories obtained with proportional navigation. 
Simulation details are as follows: the steady-state concentration field for a 
source that emits continuously with rate F = 0. 1 mM/s and diffusivity 
D= 0.05cm 2 /s is given by C(r) = F/[4iiDr). The agent starts at a random 
location Uo, yo) at r(f = 0) = 3 cm from the odor source. Its initial 
orientation is taken randomly between — jr and n. The equations of 
motion (Equation 1 as well as a = dy/dt) are integrated with a time step 
of 0.1 s. The speed of the agent is v = 0.4mm/s. The turning rate a) is 
calculated from Equation (7) with K = 3, / = 1 mm, SCi gt = Q— Q (see 
panel D) and C = (Q + C r )/2. To account for intermittent and noisy 
observations, the relative concentration change WC/ at /C) is not sampled at 
every time step but with 0.1 probability and multiplicative noise jj(<5C/ at /C) 
is added to the measurements (ij = uniform random variable with 
mean = 0 and s.d. = 0.1). 



The sensory input required for propor- 
tional navigation is examined in the next 
section. 

WHAT KIND OF SENSORY 
INFORMATION IS NEEDED FOR DIRECT 
ERROR CORRECTION AND 
PROPORTIONAL NAVIGATION? 

Drosophila larva and C. elegans assess the 
chemical gradient by comparing the stim- 
ulus intensity over time. Yet, it is not 
known exactly how the gradient is sam- 
pled. It could be sampled along the direc- 
tion of motion while the animal moves 



forward (Ferree and Lockery, 1999) or 
perpendicular to the direction of motion 
during head sweeps (Ward, 1973; lino and 
Yoshida, 2009; Gomez-Marin and Louis, 
2014). The Drosophila larva moves at 
v ~ 0.5-1 mm/s and may integrate infor- 
mation on time scales of a few sec- 
onds, and low-amplitude head sweeps lead 
to side-by-side deflections of the head 
of / ~ 0.7 mm (Gomez-Marin and Louis, 
2014). Thus, measurements along or per- 
pendicular to the direction of motion 
would be performed over distances of 
about 1 mm. Applying the same logic as 



in Gomez-Marin et al. (2010), we esti- 
mate that concentration differences are 
10-fold larger than the noise level so 
that both lateral and longitudinal sam- 
plings are viable options. In what fol- 
lows, we therefore investigate the two 
possibilities. 

In Figure 1C, we consider that the agent 
takes measurements every St sec dur- 
ing forward locomotion. The longitudinal 
gradient component is estimated along the 
direction of motion as SCi on /(vSt) where 
8Ci on is the measured difference in con- 
centration and vSt is the distance traveled 
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between two consecutive measurements. 
In the Supplementary Materials, we show 
that the longitudinal gradient component 
is related to the local bearing angle as 
follows. 

r ( 8Ci on \ 

In Figure ID, we consider side-to-side 
movements of the head sensor of ampli- 
tude I perpendicular to the direction 
of motion so that the lateral gradi- 
ent component is estimated as SQ at /l 
where 5Q flf is the difference in concen- 
tration from one side to the other. In 
the Supplementary Materials, we show 
that the lateral gradient component is 
related to the local bearing angle as 
follows. 

. r ( SC, at \ 

How to make use of this sensory informa- 
tion? From Equation (2), error correction 
requires the value of f). Taking the 
inverse cosine of Equation (5) leads 
to ±/3 whereas taking the inverse sine 
of Equation (6) leads to f} or it — /3. 
These ambiguities prevent the use of 
either the lateral or the longitudinal 
gradient component alone. In con- 
trast, an efficient implementation of 
Equation (2) using both longitudinal 
and lateral gradient components is pro- 
vided by w = K atan2(sin/3, cos/j) 
where atan2 is the four-quadrant inverse 
tangent. 

From Equation (4), proportional navi- 
gation requires the value of sin /3. On the 
one hand, the use of Equation (5) leads 
to sin p = ±y 1 — cos 2 /3 which makes 
the turning direction unknown. On the 
other hand, the lateral gradient com- 
ponent alone is sufficient for propor- 
tional navigation. From Equations (4) 
and (6), proportional navigation writes as 
follows. 

v/SC,at\ 
a) = K l{—) (7) 

We note that the turning rate a> = 
dy/dt in Equation (7) is proportional 
to the relative change in stimulus 
intensity SQ at /C so that proportional 
navigation obeys the Weber-Fechner 



law (Kandel et al, 2000). Figure IE 
shows numerical results of Equation 
(7) (simulation details in Figure cap- 
tion) that are in good agreement 
with experimental data. We also note 
that the robotic implementation of 
Equation (7) appeared to be very 
robust in real conditions (see illustra- 
tive video at http://youtu.be/XvVERq4 
h_uc). 

From experimental studies (lino and 
Yoshida, 2009; Gomez-Marin and Louis, 
2014), we know that the relationship 
between turning rate and bearing error 
in Drosophila larva and C. elegans is 
approximately sinusoidal. Also, in both 
animals, the turning rate is clearly cor- 
related with the lateral gradient compo- 
nent and not with the longitudinal gra- 
dient component. Here, we show that 
proportional navigation produces a sinu- 
soidal relationship between turning rate 
and bearing error and merely employs 
the lateral gradient component. In con- 
trast, direct error correction is linear per 
se and requires both lateral and longitudi- 
nal gradient components. Together, these 
data suggest that klinotaxis reflects pro- 
portional navigation rather than direct 
error correction. Nevertheless, because 
it is fairly impossible to rule out that 
another model may explain the data 
equally well, the conclusion that klino- 
taxis is directed through proportional nav- 
igation remains a viable hypothesis in 
the absence of compelling experimental 
evidence. 

Proportional navigation is a strategy 
employed by certain predators to track 
unpredictably moving targets, e.g., bats 
chasing insects (Ghose et al., 2006). In 
a prey pursuit, proportional navigation 
results in an advantage for the pur- 
suer that is known as motion camou- 
flage (Mizutani et al., 2003). What could 
be the interest of proportional naviga- 
tion for klinotaxis? At first sight, the 
sine relationship between turning rate 
and bearing error does not seem opti- 
mal because small turning rates are pro- 
duced for large bearing errors. However, 
the continuity of the sinus with zero cross- 
ing at ±180° may prevent chaotic reori- 
entations in the presence of noise. In 
addition, proportional navigation (unlike 
direct error correction) requires mini- 
mal sensory information; i.e., the lateral 



gradient component. Thus, Drosophila lar- 
vae and C. elegans may prefer proportional 
navigation over direct error correction 
because of its simplicity and robustness, 
just like proportional navigation is the 
preferred method of missile guidance 
because it provides the best perfor- 
mance with minimal sensory informa- 
tion (Murtaugh and Criel, 1966; Zarchan, 
2012). 
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